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Abstract 

Surface roughness is known to easily suppress the adhesion of elastic surfaces. 
Here a simple model for the contact of viscoelastic rough surfaces with significant 
levels of adhesion is presented. This approach is derived from our previous model 
[E. Barthel and G. Haiat Langmuir, 18 9362 2002] for the adhesive contact of 
viscoelastic spheres. For simplicity a simple loading/unloading history (infinitely 
fast loading and constant pull-out velocity) is assumed. The model provides 
approximate analytical expressions for the asperity response and exhibits the 
full viscoelastic adhesive contact phenomenology such as stress relaxation inside 
the contact zone and creep at the contact edges. Combining this model with 
a Greenwood- Williamson statistical modeling of rough surfaces, we propose a 
quantitative assessment of the adhesion to rough viscoelastic surfaces. We show 
that moderate viscoelasticity efficiently restores adhesion on rough surfaces over 
a wide dynamic range. 
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1 Introduction 



Surface roughness is a major parameter for the control of the adhesion of vis- 
coelastic materials and as such is central to many technological processes. As 
an example, in the automotive industry, when assembling a polymer interlayer 
with two glass sheets to form a laminated windshield, a periodic roughness is 
intentionally generated on the polymer to reduce the adhesion and facilitate po- 
sitioning before the final heat treatment [I]. Roughness also plays a role in the 
self-adhesion mechanism (or tack) of soft polymers [21 E] like repositionable ad- 
hesives. In a similar way, in the glass industry, glass molding is performed at 
elevated temperatures for easy glass flow but at the risk of adhesion. Indeed ad- 
hesion of hot glass to the mold halts production until the mold has been cleaned 
of its adherents. The roughness of the mold strongly impacts the adhesion of 
the viscoelastic glass gob [I] and a stochastic but carefully tailored roughness, 
obtained through shot peening of the steel surfaces for instance, will promote an 
easy molding process. 

The adhesive contact to rough viscoelastic surfaces is a complex problem 
however because it couples the statistical approach of the adhesive contact to 
rough surfaces with the difficulties involved in the viscoelastic adhesive contact. 
For elastic solids, the adhesive contact to rough surfaces is reasonably well under- 
stood. Building upon the Greenwood- Williamson independent asperity approach, 
Fuller and Tabor [5] have given a useful description of the strong reduction of the 
adhesive force incurred when roughness increases. They introduce an adhesion 
parameter a s /d~ c where a s is the standard deviation of the summit height dis- 
tribution and d~ c is the maximum single asperity extension before rupture. This 
adhesion parameter is a measure of the energy balance between stretched and 
compressed asperities. If the adhesion parameter is large, the compression of the 
upper asperities dominates the traction of the adhesive lower asperities which 
significantly depresses the adhesion. 

The calculations are relatively straightforward for elastic contacts which are 
essentially reversible. For viscoelastic contacts, adhesion depends upon the full 
contact history so that the calculations are more involved. In the absence of 
adhesion, the viscoelastic contact problem was first solved by Ting [BJ. Subse- 
quently, important results for the viscoelastic crack problem were obtained by 
Schapery [TJ, [S] and Greenwood and Johnson [11110]. Hui and coworkers [TTl[T2l[TB"] 
have tackled the problem of viscoelastic adhesive contacts. We recently proposed 
a model for the full viscoelastic contact of adhesive bodies [131 [T5l [16]. 

The present paper inquires into the application of this model to the adhesive 
contact to rough surfaces. To model the response of a distribution of asperities 
a simplified description of the viscoelastic contact is needed and the first part of 
the paper develops such a minimal model. It is devised so that the two main 
phenomena characteristic for viscoelastic adhesive contact - stress relaxation un- 
der the contact and creep induced by the adhesive interactions at the contact 
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edges - are preserved while keeping calculations to a minimum. The model ex- 
hibits typical viscoelastic adhesive contact phenomenology, and in particular the 
time lag between indenter retraction and the actual contact radius recession, also 
called "stick period" [H]. In the second part of the paper examples of the use 
of this approximate model for rough surfaces are presented. The impact of the 
pull-out velocity on the pull-out force for typical rough surfaces is calculated. 
The results emphasize how a moderately viscoelastic material effectively restores 
the adhesion lost through roughness. 



2.1 Single asperity contact 

Modeling the viscoelastic contact of an axisymmetric asperity in the presence of 
adhesion is complex because, as in all viscoelastic problems, the solution depends 
upon the full contact history. 

In addition, for contact problems, one must deal with mixed boundary con- 
ditions: the surface normal displacement u(r) is prescribed inside the contact 
zone whereas the normal surface stress a(r) due to the adhesive interactions is 
specified inside the so-called cohesive zone (Figure [T]), outside the contact zone. 
Moreover the boundary between these two regions is not fixed but moves during 
the contact history. For a viscoelastic material, the solution is reasonably simple 
for growing contact radius [HI [TTJ . However, during pull-out (i.e. decreasing 
contact area), the residual deformation has to be taken into account. This in 
turn depends upon the full stress history, which leads to a complex situation, 
both for adhesionless [6] and adhesive [HJ US, [13j HH [15] contacts. 

2.2 Resolution method 

Under the assumption of axial symmetry, we resort to two auxiliary functions g(r) 
and 6(r), which are the following Abel transforms of respectively the distributions 
of normal surface stress a(r) and normal surface displacement u{r): 



2.2.1 Boundary conditions 

A first benefit of these transforms is their pertinence for axisymmetric boundary 
conditions. In particular, inside the contact zone (r < a), the normal surface 
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displacement is known through the contact condition 

u(r) = 5- h(r) (3) 

where 5 is the penetration and h(r) the shape of the indenter (Figured]). Then, 
using Eq. [21 inside the contact zone 

6(r) = 5-5 (r) (4) 

where 

Mr)= «r ds ^L 2 . (5 ) 

dr Jo Vr 2 - s 2 

It turns out that So{r) is the Hertz penetration for a contact radius r as shown 
in the next section. Note that this function So{r) is determined by the shape of 
the indenter only. It is equal to r 2 /R for a sphere of radius R and to |r/ tan/3 
for a cone of apical angle (3. 

Conversely, as mentioned above, the adhesive stress distribution <r(r) outside 
the contact zone is known. More precisely, in the present case, we will show in 
section 12.2.31 that it can be self-consistently determined through an independent 
set of equations, under the assumption of small cohesive zone size. Therefore the 
g function is known outside the contact zone (r > a). 



2.2.2 Mechanical Equilibrium — Viscoelastic materials 

The second benefit of these auxiliary functions is that they naturally handle the 
long range nature of the elastic field: for instance, under the conditions of linear 
elastic behaviour (Young's modulus E, Poisson ratio v) and absence of shear 
stresses at the contact, mechanical equilibrium leads to 

9(r) = fV) (6) 



where 



The contact problem is solved under the assumption of continuity of the stress 
distribution at a. The penetration is then directly obtained by 



g(a) = *L-0(a) (8) 

For example for an adhesionless contact, cr(r) = for r > a so that g(a) = 
and 5 = 5o(a) as mentioned in section 12.2.11 If adhesion is present, g(a) is not 
zero and a more complex situation arises. 
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For a linear viscoelastic material, under the approximation of constant Poisson 
ratio [IB] we introduce the creep and the relaxation functions 

<P(t) = ^-Jit) (9) 

m = (io) 

with 

fct) = l + (11) 

4>(t) = k + (l-k)e% (12) 

where /i is the creep time. The parameter k lies between 1 (elastic) and 
(Maxwell). Due to the decoupling between spatial and temporal responses, the 
equilibrium equation (Eq. [6]) will now be 

g(r,t)= fdnl){t-T)-^0{r,T) (13) 

JO CLT 

or its inverse 

0(r,t) = J*dT(f>(t-T)-^-g(r,T) (14) 

where the relaxation function tp(t) and the creep function <f)(t) are inverse for this 
product of convolution. In the case of increasing contact radius, Eqs. H] and HH 
determine the penetration by 

6(t) = S (a(t)) + f dr<f>{t - r)^-g(a(t), r) (15) 
Jo or 

This result is equivalent to previous formulations [T2l [T7] . 

The decreasing contact radius case is more intricate. Taking into account the 
domains where g and 9 are known respectively, we have to resort to Eqs. HI [T3] 
and HH to determine the penetration by the integral equation 

g(a(t),t)= f dr^(t-r)^-{6(r)-S (a(t))} + g(a(t),t) (16) 

Here t a -(a(t)) is the time at which the contact radius was first equal to a(t), 
during the loading phase. The corrective term 

rt a -(a(t)) d ( r , , d ,\ 

g(a(t),t) = j o dr^t-r)— (jf dr> >(r - r')—g{a{t),r')^ (17) 

is unessential, except for the cross-over between increasing and decreasing contact 
radius regimes (a = a max , t = t max ). 
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Then for a given history of the penetration, one can calculate the history of 
the contact radius and subsequently the force [TUESHS]- Under the assumption 
of small cohesive zone size (see section 12.2.31 for the cohesive zone and section [6] 
for details of the calculation) the force is given by Eq. 19 in [1] 



d /■min(a(t),a(r)) 

P(t) = 4 drip{t - t)— dr(5(T)-5 (r)) (18) 



o dr 

Note that this equation is valid for both inward and outward runs. Eqs. [T6land[T8l 
involve stress relaxation inside the contact zone but also crack tip viscoelasticity 
through g(a(t),t), as discussed in more details in the next paragraph. 

2.2.3 Cohesive Zone 

For any adhesive contact, in the context of the present (Sneddon) method, the 
quantity g(a), which is the function g evaluated at the edge of the contact zone, 
plays a central role. For the contact zone, it determines Eqs. H51 [TH] and [TH] 
between macroscopic contact variables. But we have previously shown that, for 
the cohesive zone, g(a) determines the stress intensity factor [19] characteristic 
of the local deformation process at the crack tip. 

Here we dicuss, for the specific case of a viscoelastic adhesive contact, how 
g(a(t),t) is determined by the cohesive zone parameters. This results in an inde- 
pendent relation between g(a) and the contact radius velocity a through Eqs. |2"U| 
EH E2 and El 

This problem is similar to a viscoelastic crack [7J EJ EJ [10] and calls for a 
consistent treatment of an adhesive-stress induced creep inside a cohesive zone of 
finite size (c(£) > r > a(t), Figure [JJ). The equivalence of various methods for the 
description of the cohesive zone has been demonstrated earlier [20] • In brief, the 
details of the description of the stresses inside the cohesive zone is unimportant 
as long as the effective range of these interactions is preserved. Greenwood and 
Johnson [21] have proposed a description (the so-called " Double- Hertz" ) which is 
analytically simple and which we have retained for our viscoelastic models |14j . 
Following these authors, we assume a specific stress distribution which results in 

71- 7< 2 — 

«c> - r°w=r- (19) 

Assuming that the size of the cohesive zone is small compared to the contact 
zone (a 3> c — a = e) and taking the limit when r —> a gives: 

g(a) = -^a V2a~e (20) 

For a decreasing contact radius the self-consistent equation coupling the total 
adhesion energy and the deformation of the contact zone writes (Eq. 23 of Ref. 

!): 

w = ^a 2 eMtr) (21) 
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where w is the adhesion energy, 

= ! J^drr^t-r) (22) 

and t r is a fracture dwell time. For non vanishing contact radius velocity, this is 
the time it takes for the crack tip to move a distance equal to the cohesive zone 
size. It is given by: 

t r (t) = =^ (23) 
a[t) 

The simultaneous resolution of Eqs. l2"UH2"3l determines the velocity dependence 
of g(a). It is typical of viscoelastic crack problems : compare our set of Eqs. |2"U| 
[2T1 and [231 with, for instance, Eqs. 1, 47 and 48 in [7j. As a result, g(a) increases 
with d as displayed on Figure EJ again compare e. g. with Fig. 6 b in [TO] . 



2.2.4 General Model stick period 

For the viscoelastic adhesive contact, the main phenomenon we have identified 
[HI [15] is the "pinning" of the contact upon unloading (Figure [3]). More precisely, 
for such viscoelastic contacts, there is a typical time lag between the retraction 
of the indenter and the actual contact radius recession. We have called "stick 
period" this period over which the contact radius stays pinned, before contact 
edge retraction sets in, finally leading to contact rupture. 

The physical interpretation of the stick period phenomenon has already been 
given in |14j . Relaxation of the stresses inside the contact zone reduces the stress 
intensity factor and hence g(a). This results in a slow crack motion (Figure [2]) 
and the contact radius practically comes to a halt. When pulling the indenter 
out, enough tensile stresses must be re-built inside the contact zone before g(a) 
is large enough for the contact edge to be restored into propagating condition 
(Figure EJ), and this depends upon the competition between stress relaxation and 
pull out velocity as discussed in [T6j . 



2.3 Approximate Model 

We need a simple description of the viscoelastic adhesive contact for rough surface 
calculations. We restrict the developments to the spherical asperity (radius R) 
usual in rough surface models although extension to other asperity shapes should 
be straightformard. We also introduce two major simplifications to the full model. 
One is an assumption on the type of loading/unloading. The other one is an 
approximation for the stick period. 



2.3.1 Simplifying assumption — Direct unloading 

We consider fixed-grip loading, i.e. a loading where the penetration is prescribed. 
In addition, the loading phase is considered to be infinitely fast and the pull-out 
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velocity 5 a (negative) constant (Figure ED, so that the penetration history simply 
writes: 

5(t) = (5 + t5)Q(t) (24) 

where 5q is the initial penetration^, taken positive when the two surfaces are 
pushed one into another. is the Heaviside step function. Due to the infinitely 
fast loading phase the quantity t a ^(a(t)) defined in the unloading phase (sec- 
tion I2.2.2H is equal to zero. 



2.3.2 Stick period approximation 

This approximation stems from the observation that for large penetrations, the 
contact radius upon loading is primarily determined by the instantaneous elastic 
response and that crack tip creep during loading is not a dominant process |14j . 

As a result, in the following, the stick period is approximately modeled by 
a period of constant contact radius prescribed by the instantaneous modulus. 
This will be a good approximation except for small initial contact radii where 
the adhesive term is dominant. 

In summary the approximate model is as follows: 

1. calculate the initial contact radius from an elastic model with the instan- 
taneous Young's modulus; assume this contact radius stays constant to the 
end of the stick period. The main viscoelastic process at work here is the 
stress relaxation inside the contact zone, which directly impacts the force 
(section |2~3~3"1) : 

2. determine if there is a stick period and if so calculate ending time tf (sec- 
tion E33D; 

3. calculate the unloading phase (section [2. 3. 5p and in particular the adhesion 
force. 



2.3.3 Stick Period - Force 

The force in the stick period is obtained from Eq. [18] taking into account the 
discontinuity of the penetration at r = 0. This leads to (section 16.11) : 

Pi(t) = m) (ff - g °^f* ) + [ dr^t - r) (25) 

which can be rewritten as: 

Pi(t) = $(t)P„ + 4a 5 f dnl){t - r). (26) 

Jo 

1 This initial penetration cannot be confused with the 6o(r) function introduced earlier and 
which will now be substituted by its specific value for a sphere r 2 / R in the rest of the paper. 
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where the normalized form of ip (Eq. [T2l has been used. P is the initial force 
and 6q the cohesive zone size at t — which can be calculated from Eqs. |2"U1I2"T1 
The first term results from the relaxation of the initial force and the second one is 
due to the flat punch viscoelastic response and is proportional to the (constant) 
pull-out velocity. 



2.3.4 Approximate determination of the end of the stick period 

At the end of the stick period a{tj) = a because the contact radius increase 
has been neglected by assumption. However we do take into account the crack 
tip evolution which is reflected in a cohesive zone size e(tf) smaller than e . 
Since the contact radius is a constant, the crack dwell time t r (section 12.2.31) is 
approximated as the crap tip age, which is also the contact age. Therefore, when 
the contact radius starts decreasing at t — tf, the dwell time t r is taken equal 
to tf. Inserting Eq. [2H into Eq. |2"U1 provides an approximate criterion for the 
begining of the contact radius decrease: 



8 « ttM= im (27) 

When the contact radius decreases, Eq. [161 becomes 

g(a(t),t) = 6 jf*dr0(t - r) + ^(t) (<5 - ^pj (28) 
Then tf is solution to the following implicit equation: 




dTi'(t, - t) + ./-^ = (29) 



Note that the purely elastic case is recovered for t = 0. It is also recovered 
for constant ip{t) (elastic material) because 5(t) = 5o + 5t, and for long waiting 
times where ip(t) ~ ^(+oo) and Jq^^cZt ~ t^(+oo). 

Injecting the initial elastic solution, tf is identified as the non zero root of 
function /, with 

f(t) = 6 dr^(t - r) + v^F^ ^ -j== - m\[-L) (30) 



Time t / only depends on the mechanical properties of the material and on the 
history of the penetration (So and S). Considering a given contact experiment on 
a rough surface, the mechanical and adhesive properties as well as the pull-out 
velocity are identical for all asperities. Therefore, tf will only be a function of the 
initial penetration <5 which is depends on the initial height of a given asperity. 
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The numerical calculation of tf(8o) is achieved using a bisection method. From 
Eq. [301 and knowing that a increases with 8q, t/(5o) is a strictly increasing 
function of 5o and can thus be numerically inverted to obtain the function Zf(t), 
representing the initial penetration necessary for a given asperity to reach the 
end of the stick period at time t. This function will be useful for the computation 
of the total force contribution (section I2.5P from the asperity height distribution. 

2.3.5 Decreasing contact radius 

The computation of the load i n the case of a decreasing contact radius is 

performed using Eq. [TH] which leads (section 16.21) to: 

P 2 (t) =4a(t)m {^o-^pj +4a(t)6£dTil>(t-T) (31) 

The two terms are respectively the elastic force for a spherical asperity with 
the relaxed modulus at time t in presence of adhesion and the response of the 
viscoelastic medium to a flat punch displacement (contact radius a(t)). Note that 
the purely elastic case is recovered for elastic materials because S(t) = 5 + St. 
Eqs. EU [23J EH] and EH completely describe the decreasing contact radius phase. 



2.4 Single Asperity — Normalization 

For a spherical asperity, following Maugis [M1I22], the normalization of the con- 
tact variables is performed by introducing the following quantities: 

P = (32) 
A = ~, TTTTv (33) 



( 3nwR 2 \ 


1/3' 


\ 4E* J 




5 




( 9tt 2 w 2 R 


j 1/3' 






2a 





TTWWE* 2 ^/ 3 



9R 



(34) 
(35) 



G(A) = -^\V2M. (36) 

were e is the normalized size of the cohesive zone (same length scale as contact 
radius a). In normalized form, g(a) is given by: 



7T 



G(A ) = --Xy/2A e (37) 
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In normalized form, the self-consistent equation is given by: 

1 = ^\ 2 ~e(t)Mtr(t)) (31 

where <f>\ is given by: 

2 r* 



MU) = -2 / rfr(r)0(r) (39) 
t r Jo 

Consequently, Eq. [2H1 used for the determination of the end of the first regime 
now writes: 



m(m-4}+5 d<nJ>(t f -r) + 2 , -J7— = (40) 

In the stick period, the force for one isolated asperity (Eq. [TBI) is now given by: 
AW = m (A 3 (t) - A l 5 V0) + \a,~5 £ <*r^(* - r) (41) 



2 

When the contact radius starts decreasing, G(A(t)) (Eq. |2"81) writes as follow 



G(A(t)) = 5 y dr${t - r) + ^(i) (5(0) - A 2 (t)) (42) 
and the force in the second regime is expressed through: 

h(t) = m (~A(t)~5 - + ?A(t) S jf - r) (43) 

2.5 Accounting for a rough surface 

We follow the roughness description proposed by Greenwood and Willamson [23] 
in their adhesionless elastic contact study. Namely, a surface of nominal area Aq 
with N asperities is in contact with a rigid flat plane. All asperities are assumed 
to have the same radius of curvature R and a Gaussian height distribution of 
standard deviation <t«: 



's 

.2 



X(z) = 7= exp I (44) 



a s V2n V 20 s, 

where x( z ) is the probability that the summit of an asperity stands between z and 
z + dz. For a penetration of the rigid flat 5q , the number of asperities in contact 
is given by: n = N dz\(z) and all contact variables can be obtained with the 
same integration. Similarly, the calculation of the load at time t is performed by 
integrating the viscoelastic response of each asperity initially in contact: 

P(t)=N dz X (z + 5o)P a (z,t) (45) 
Jo 

where P a (z, t) is the load at time t of an asperity undergoing an initial penetration 
of z — 5 . 



12 



3 Results and discussion 



3.1 Single asperity — Impact of pull-out velocity 

Force vs penetration curves (force plots) for a single asperity were calculated with 
A=5 - a value consistent with our assumption of small but finite cohesive zone 
size - and for a moderately viscoelastic material with k = 0.5. Figures H] and [5] 

display results obtained for So = 5 for pull-out velocities S ranging from —0.01 to 
-50. 

For very low pull-out velocity values (\5\ < 0.1), there is an initial rapid 
decrease of the force, followed by a linear variation of the force as a function of 
penetration. This is due to the fast relaxation of the (predominantly compressive) 
stresses. The stress field is then controlled by the relaxed modulus which leads 
to constant contact stiffness and a linear force- distance curve segment, before 
the contact radius starts decreasing. When the penetration decreases further, 
the force becomes tensile (negative), displays an incurvation, levels off, finally 
increases and jumps to zero: this evolution of the force results from contact 
edge recession and asperity snap-off. The absolute value of the maximum tensile 
force defines the pull-out force F pu u out , and the (positive) effective adhesion w e ff 
through w eff = F pullout /(3/2 ttR). 

The very low velocity case is an elastic adhesive contact with an effective 
modulus equal to the relaxed modulus and a pull-out force equal to 3/2 which 
is the JKR value. At slightly larger pull-out velocities, we observe that the pull 
out force increases with pull-out velocity from this 3/2 value: higher pull-out ve- 
locity leads to increased contact radius velocity. Enhancement of the viscoelastic 
crack tip dissipation and effective adhesion follow as amply demonstrated in the 
literature [7J El El [HI [121 El El] . It results from the combination of a fully relaxed 
contact zone and a less relaxed cohesive zone. 

In contrast, for large pull-out velocities (\5\ > 1), the initial precipitous drop 
of the force is no longer present and the pull-out force decreases with pull-out 

velocity. Eventually for very large pull-out velocity (e.g. 5 = —50) the behavior 
tends toward an elastic behavior and the pull out force reverts to 3/2. Indeed, for 
large pull-out velocities, the relaxation of the compressive stresses in the contact 
zone is far from complete during the experimental time. The pull-out force is 
lower than for low pull-out velocity because the compressive contact stresses 
favor contact rupture. 

As a consequence, there exists a transition with an optimum pull-out velocity 

(here of around 5 = — 1) for which the pull-out force is maximum. This result 
specifically originates from the consistent inclusion of both viscoelastic crack tip 
dissipation and contact stress relaxation. 
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3.2 Rough surfaces 



As an example of results for rough surfaces, the same material parameters A = 5 
and k = 0.5 are chosen. The initial penetration is 5q = 5. 

3.2.1 Force distribution 

Normalizing the roughness as the penetration (Eq. IH4"|) . we compute the force 
curve for a rough surface with a normalized standard deviation d s = 2. Com- 
puting the contributions of asperities displaying each regime (stick period vs. 

decreasing contact radius) to the total force for pull-out velocity 5 = — 1 (Figure 
we observe that the proportion of asperities initially in contact is almost 100%, 
as expected for such a ratio of initial penetration 5 to roughness <f s . This type of 
contact is relevant to glass molding where a rather homogeneous distribution of 
asperities, with narrow asperity height and radius distributions, is generated on 
the mold surface. Note that the Greenwood- Williamson model is particularly apt 
at treating such a surface morphology. It would be less adequate for a less regular 
surface with a wide distribution of asperity height and/or a wide distribution of 
asperity radius. 

The force is initially dominated by asperities pinned in the stick period (black 
dashed line) . A gradual transition where the contact zones of the lower asperities 
de-pin and start to recede is observed: the contribution of pinned asperities 
drops to zero. Then the lowest asperities start to undergo contact rupture. It is 
in the middle of this transition that the maximum tensile force (pull-out force) 
is recorded. 

3.2.2 Impact of roughness — Constant pull-out velocity 

Force plots with increasing roughness (from d s = 0.1 to <f s = 5) are displayed 
on Figure [7J For small d s , the results track the single asperity limit: all the 
asperities in the population exhibit an identical behaviour, having nearly the same 
initial penetration. For larger a s , the initial contact force increases as expected 
from the elastic (JKR) limit: for a single asperity the increase of the force with 
initial penetration is more rapid than linear. Simultaneously, the duration of the 
contact increases and the pull-out force decreases with roughness as in the elastic 
case [5]. Indeed, for one isolated asperity, the pull-out time increases with initial 
penetration. Therefore, when a s increases, the distribution of the pull-out times 
spreads out and the different asperities contribute less constructively to the total 
pull-out force, which therefore decreases. 

3.2.3 Impact of pull-out velocity — Constant roughness 

Figures El and [9] display force plots at different pull-out velocity values, from 
5 = —0.01 to 5 = —10 for a normalized standard deviation d s = 2. The overall 
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picture is similar to the single asperity case (section l3.ip and the same transition 
from low to high velocity regimes is observed. At low pull-out velocities, a fast 
decrease of the force due to stress relaxation is followed by a quasi elastic behavior 
while at large velocities, the system is purely elastic. The pull-out forces in these 
two limit regimes are differentiated by the roughness. This is because the elastic 
moduli involved are different: indeed, the relaxed modulus is pertinent for the 
very low pull-out velocity regime, and the instantaneous modulus for the very 
high velocity regime. A larger elastic modulus will result in a larger adhesion 
parameter and a lower adhesion force [5], [25] . A pull-out force enhancement is 
recorded between these two limit regimes because of the combination of partial 
contact stress relaxation and crack tip dissipation. 

It is interesting to assess the details of the transition between these two 
regimes. Calculation of pull-out forces as a function of pull-out velocities were 
performed for a rough surface with a roughness of a s = 2, a cohesive zone pa- 
rameter A = 5 and a moderately viscoelastic material k = 0.5. The pull-out force 
as a function of pull-out velocity has been normalized to the rough surface limit 
in the high velocity elastic regime (instantaneous modulus). The results are dis- 
played on Fig. [TU1 where the isolated asperity result is also shown for comparison. 
As expected both exhibit enhancement of the pull-out force in the intermedi- 
ate regime where contact stress relaxation and crack tip creep occur at the same 
time. Figure [TU] clearly evidences that the enhancement of adhesion is much more 
pronounced on rough viscoelastic surfaces than for an isolated asperity. Indeed 
for elastic materials a roughness as large as a s = 2 severely depresses the pull-out 
force due to the stored elastic energy of the higher asperities. However for a 
viscoelastic material these elastic stresses rapidly relax while the single asperity 
adhesive force is enhanced through crack tip dissipation. As a result, even a very 
moderately viscoelastic material is sufficient to enhance the effective adhesion on 
rough surfaces by an order of magnitude over a large range of pull out forces 
(Figure [TO]) . 

4 Conclusion 

A simple model for the adhesive contact of rough viscoelastic surfaces has been 
developed based on a Greenwood- Williamson roughness distribution. For that 
purpose a simplified model has been proposed for the adhesive contact of vis- 
coelastic asperities, assuming instantaneous loading and constant pull-out veloc- 
ity. It retains the full phenomenology of the viscoelastic contact including stress 
relaxation inside the contact zone and creep in the cohesive zone. 

For elastic surfaces, it is well known that increasing roughness rapidly sup- 
presses adhesion. Our calculations give a quantitative assessment of how effi- 
ciently viscoelasticity restores adhesion on rough surfaces through the interplay 
between contact stress relaxation and cohesive zone creep. The phenomenon 
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operates in a wide range of pull-out velocities. 
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6 Appendix — Computation of the total force 

The total force writes as [15] : 

r+oo 

P(t) = 4 / g(r,t)dr (46) 
Jo 

The contribution of the cohesive zone is of higher order in e (cohesive zone size) 
and is neglected, which is typical for large A. Then Eq. [46] becomes P(t) = 
4 Jq^ g(r, t)dr and Eq. [TBI becomes 

g(r 1 t)=[ t d T^(t-T)^9(r,r) (47) 
Jt a -(a(t)) dr 

with #(r, r) = 5(t) — 5o(r) which is Eq. HJ Inserting Eq. 1471 into 1461 and inverting 
the order of integration (cf also Ref. [15], section 3.2.4) results in Eq. [T8] 



6.1 During the stick period 

During the stick period by assumption a(t) = ao is constant and min(a(t), a(r)) = 
ao- In addition for a sphere of radius R loaded as specified in section 12.2.11 
combining Eqs. H] and [21] 

e( r , r) = (5 + t5- r 2 /R)0{r) (48) 

Then Eq. [TS] gives: 

Pi(t) = 4^(t) / dr(S - — ) + 4a 5 / dr^{t - r) (49) 
Jo R Jo 

P 1 (t) = 4^(t)Wo-^)+4a o( 5 fdr^t-r) (50) 

SR Jo 

where Pi denotes the force during the stick period. The initial penetration Sq is 
given by the double Hertz model and writes: 



a ao-KJa eo 

"o = t= • (51) 

R V2E* v ; 

This last equation, in combination with the previous one, gives Eq. 
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6.2 Decreasing contact radius 

By assumption the contact radius is constant or monotonously decreasing from 
r = so that min(a(t), a(r)) = a(t) and Eq. [181 directly results in Eq. [3TJ 
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Captions and Figures 

Fig. [D Schematics of the viscoelastic adhesive contact showing the gap between 
the surfaces and the normal surface stress distribution. Special emphasis is given 
to the cohesive zone where the attractive interactions operate across the gap 
between surfaces. 

Fig. [2) Schematics of the relation between g(a) - proportional to the stress 
intensity factor - and the crack velocity (after [T4"]). 

Fig. [3j Prescribed penetration as a function of time {5(t)) and resulting con- 
tact radius history (a(t)). The constant radius approximation proposed in the 
present model is also shown as a appro {t). The end of the stick period is t/. 

Fig. |U Single asperity force plot for different values of the pull-out velocity 

(from 5 = —0.01 to 5 = —1) and for the elastic limit (JKR) with the relaxed 
modulus. 

Fig. EJ Single asperity force plot for different values of the pull-out velocity 

(from 5 = — 1 to 5 = —50) and for the elastic limit (JKR) with the instantaneous 
modulus. 

Fig. El Force distribution during adhesion rupture on a rough surface showing 
the contributions from asperities in the stick period and asperities with receding 
contact radius. The proportion of asperities in contact is shown on the right hand 
axis. 

Fig. [Tj Force plots for rough surfaces with increasing roughness. Here Sq = 5 
and 6 = — 1. 

Fig. [HJ Force plots for a rough surface with o~ s = 2 for increasing pull-out 
velocity at low velocities. Also shown is the elastic limit for the relaxed modulus. 

Fig. [9j Force plots for a rough surface with o~ s = 2 for increasing pull-out 
velocity at high velocities. Also shown is the elastic limit for the instantaneous 
modulus. 

Fig. [101 Pull-out force enhancement over elastic function of pull-out 

velocity for an isolated asperity and for a rough surface with a s = 2. 
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